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We will present here the simulations of RICH detectors which will be used in the CKM experiment. We will verify their 
performance, critical to the experiment. 
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En este trabajo, se presenta una simulacion con el fin de verificar el desempeno de dos detectores RICH que se usaran en el 
experimento CKM. 
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1. Introduction 

In this article we present results on simulating a RICH 
detector which will be used within the CKM experi- 
ment [1], a rare kaon decays experiment recently ap- 
proved at Fermilab. The main objective is to verify if 
the resolution of the current design of the detector is 
sufficient for the CKM experiment. 

After a short description of the CKM experiment, we 
will detail the requirements for the detector system, de- 
scribe the current design, our method of simulation, and 
the results obtained so far [2]. 

2. Description of the CKM Experiment 

CKM stands for "Charged Kaons at the Main Injec- 
tor" . The main goal of this experiment is to measure the 
branching ratio of K + — > ir + vv to about 10 % statistical 
precision. This measurement can be used to extract in 
a theoretically clean way the magnitude of Vtd, an ele- 
ment of the Cabibbo Kobayashi Maskawa matrix, with 
an overall precision of 10 % [3] . 

This measurement plays an important role in testing 
the description of CV violation in the Standard Model, 
and is complementary to the current efforts of the B- 
factories 

The prediction of the Standard model [3] for the 
branching ratio of K + — > ir + vv is B(K + — > tt + vv) = 
0.75±0.29 x 10~ 10 and was first calculated by Inami and 
Lim [4]. 

The E787 experiment at Brookhaven, with a stopped 
Kaon beam, recently observed two clean events [5], lead- 
ing to a measured branching ratio of B = 1. 57^0 82 ' 
10 -10 . A follow-up experiment at Brookhaven, E949, 
will exploit fully the stopped kaon beam technique, and 
expects to observe 5-10 events. 



To reach the 100 event level, CKM will use kaon de- 
cays in flight. Due to kinematic backgrounds, the ac- 
ceptance will only be in the 1 — 2 % range, and this, 
together with the spill structure of the Main Injector 
(1 s out of 3 s) and a running time of 2 years, requires at 
least 30,000,000 kaons passing the detector per second. 
This places high performance criteria on the detectors. 
With an expected branching ratio of 10~ 10 and 100 sig- 
nal events with less than 10 background events to ob- 
serve, the experiment has the challenge of suppressing 
all sources of possible backgrounds to the 10 -12 level. 

The most obvious backgrounds are the main K + de- 
cays modes: K + — > n + v^ (B = 64%) and K + — > 
7T + 7r° (B = 21%). In the first, miss-identifying the fi + 
as 7r + , and in the second not registering the two pho- 
tons from the w° decay lead to identical topologies like 
the signal K + —* tt + vv. 

We will describe here the experimental method on the 
7T + 7r° mode, other modes work the same. With a very 
efficient Photo- Veto-System [9] CKM expects to reject 
this mode by 10 -7 . For the additional rejection down 
to the 10~ 12 level we use kinematic rejection: By mea- 
suring the invariant mass of the system the ir + recoils 
(the missing mass), we obtain a peak at the ir° mass in 
the 7T + 7r° two-body decay, and a continuous distribution 
in the signal mode (a three-body decay). By exclud- 
ing the mass region around the 7r°, we obtain our signal 
events. The missing mass M m i ss is given (in a very good 
approximation) by 

M miss = "^(i -P-k/Pk) +ml(l-p K /pi r ) -PkPkO 2 

where vhk is the mass of the K + , m n is the mass of the 
7r + , pk and p n are their corresponding momenta, and 
is the decay angle. These are the only experimentally 
accessible quantities in the signal decay. 

CKM will use two different types of spectrometer 
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systems to measure the momentum vectors of the K + 
and the 7r + , a magnetic spectrometer with high-rate 
MPWCs and straw-tubes, and a velocity spectrome- 
ter [6], a novel concept, consisting of two Ring Imaging 
Cherenkov counters [7]. The resolutions in momentum 
and angle of these detectors has to be good enough to 
allow for a narrow enough missing mass peak of the 7r° 
not to overlap too much with the signal region. It also 
requires a magnetic field homogeneity of better than a 
few per mil. 



3. Requirements and Design Criteria for 
the Velocity Spectrometer 

The velocity spectrometer has to measure independently 
of the conventional magnetic spectrometer the vector ve- 
locity of the incoming kaon and the outgoing pion. In a 
first approximation for designing the experiment [8] we 
assumed momentum resolutions of 0.5% and 1% respec- 
tively. 

The RICH for measuring the outgoing 7r + was mod- 
eled after the very successful phototube based SELEX 
RICH detector [10], keeping the diameter of the vessel 
about the same, but doubling the vessel length to 20 m 
and correspondingly the radius of curvature of the spher- 
ical mirror to 40 m. As radiator gas we will use Neon, 
as in the SELEX RICH. The phototubes are sensitive 
down to 160 nm, requiring also sufficient reflectivity of 
the mirrors in the VUV. 

The Kaon RICH, to measure the velocity vector of 
the incoming K + , and also to discriminate against other 
beam particles (mostly 7r + ), had to be designed differ- 
ently. We will still keep the main feature from the SE- 
LEX RICH, the use of photomultipliers for the detection 
of the Cherenkov photons. 

The beam in CKM is very parallel, and has a diameter 
of only 10 cm. Since the detector is located just before 
the decay volume, we have to avoid interactions in de- 
tector material, asking for as less as possible amount of 
material in the beam. To reduce the effect of the pho- 
totube size on the overall resolution, we also decided 
to double the focal length, leading to a doubled ring ra- 
dius. All this together lead us to a design shown in fig. 1. 
The vessel has a length of 10 m, but we double- fold the 
light-path, with a (thin) plane mirror in the beam, a 
spherical mirror of 40 m radius (20 m focal length) and 
another plane mirror, both outside the beam. 

A detailed design of the mirror arrangement is shown 
in fig. 2. The two plane mirrors have the same distance 
from the spherical mirror, as well horizontally as verti- 
cally, so that the tilt-angles of the two plan mirrors have 
the same absolute value. This assures normal incidence 
on the photo-tube plane without tilting it. With simple 
geometry it can be seen that the inclination angle a for 
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Figure 1: Principle design of the Kaon RICH for CKM 
used in this study 
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Figure 2: Detailed layout of the mirror geometry for the 
Kaon RICH. 

the plane mirrors is given by 

1 1 fA\ 
a = —6 = — arctan I — I 

2 2 V / 

where the focal length of the spherical mirror is given by 
F = L + L/ cos9. 

The size of the mirrors can be determined from the 
Cherenkov angle, given by [11] 

(1) cos 

where (3 — v/c is the velocity of the particle as a fraction 
of the velocity of light in vacuum, and n is the refrac- 
tive index of the medium. The sizes of the 3 mirrors 
are determined assuming a K + with a momentum of 
22 GeV/c, and the size of the photocathode (a matrix of 
1/2 inch photomultipliers) is determined assuming a 7r + 
with the same momentum. The sizes obtained depend 
on the refractive index n of the radiator gas. In table I 
we show typical values for the final parameters obtained 
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Diameter 


Thin plane mirror 
Second plane mirror 
Spherical mirror 
Phototube array 


43.5 cm 
90.5 cm 
79.5 cm 
90.5 cm 



Table I: Sizes for the mirrors and the phototube array 
for CF 4 at 0.8 atm. 




Pressure [atm] Pressure [atm] Pressure [atm] 



in this study. 

4. Data Analysis and Results 

After coding the full geometrical description of all rel- 
evant materials within the GEANT [12] package, and 
incorporating all the routines into the standard exper- 
iment Monte Carlo program, we studied as a first ap- 
plication the momentum resolution of the detector as 
a function of pressure for two different radiator gases, 
namely N 2 and CF 4 . N 2 was selected because it's re- 
fractive index at 1 atm of n — 1 s=s 300 • 10~ 6 gives the 
correct threshold for kaons of 22 GeV/c, but we were 
afraid of chromatic dispersion (change of n as a func- 
tion of wavelength [13]). For this reason we selected also 
CF4, with lower dispersion [14], at the cost of running 
at a non-atmospheric pressure. 

Using a 22 GcV K+, where GEANT takes care of all 
possible interactions like multiple scattering etc., includ- 
ing the production and tracking of Cherenkov photons 
and simulation of the wavelength depending efficiency 
function of the phototubes [15], we obtained a list of 
hit phototubes for every K + we simulated. For every 
event, we performed a simple ring-fit [16] to determine 
the radius of the Cherenkov ring, obtaining at the end 
an average radius r and its standard deviation ay. The 
radius r in a RICH as function of the momentum p for 
a particle with mass m is given by (with a small-angle 
approximation of eq. 1) 



(2) ^'M'-mT 

and its derivative by 

(1 - m 2 /p 2 )~ i - 



m 2 /p 2 ) 2 



Figure 3: Example results for values obtained from 
the simulation with CF 4 as radiator gas, for K + with 
22 GeV/c, all as a function of pressure. Left: Cherenkov 
Ring Radius. Middle: Ring Radius Resolution oy. 
Right: Momentum Resolution a p . 



In fig. 4 we show the final result of this study, the 
momentum resolution a p as a function of pressure for 




Figure 4: Results obtainted with the GEANT simulation 
for the momentum resolution u p for the Kaon RICH, 
with (left) CF 4 and (right) N 2 as radiator gas. The 
arrows indicate the proposed momentum resolution of 
0.5% for a 22GeV/c K+ . 



two possible radiator gases, CF4 and N 2 . The arrow 
indicates the required momentum resolution of 0.5 % for 
a 22GeV/c K + . As seen, the required resolution can 
be obtained by using CF4 as radiator gas. As shown in 
table II, we analysed the different contributions to the 
resolution. The main factor is the chromatic dispersion. 



R 
2\/2 



1 - 



1-^ (m 2 /p 3 ) 



dr 
dp 
(3) 

where R denotes the radius of curvature of the spherical 
mirror. With the simple relation a p = ^o> we can 
determine the momentum resolution a p for any set of 
parameters we wish to simulate. 

As an example we show in fig. 3 the radius r, the 
radius resolution ay, and the momentum resolution a p 
as a function of the pressure for CF 4 as radiator gas. 
An increase in pressure corresponds to an increase in 
the refractive index n of the radiator gas. 



Type 


ay 


Multiple Scattering 
Chromatic Dispersion 


0.09 cm 
0.49 cm 


Total 


0.49 cm 



Table II: Contributions to the ring radius resolution of 
the Kaon RICH with CF4 at 0.65 atm as radiator gas. 



5. Conclusions 

In this work we performed a GEANT simulation of a 
RICH detector to determine its response and momen- 
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turn resolution as a function of pressure for two different 
radiator gases. 

We obtained as a result that with CF4 at 0.8 atm pres- 
sure as radiator gas we can achieve the required resolu- 
tion of 0.5%. 

The programs and routines written for this study were 
used extensively for other design studies for the second 
edition of the proposal for the CKM experiment [9] . The 
experiment was approved by the Fermilab Director in 
June 2001. We are now improving and testing the design 
of the RICH detectors in preparation for a Technical 
Design Report. 
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